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We review an "exact semiclassical" resolution metliod for the general 
stationary ID Schrodinger equation 

- ^V'(g) + [V{q) + A]^(g) = (g G M) (1) 

with a polynomial potential 

V{q) = + Viq^~^ H h VN_iq, denoting (wi, ■ ■ ■ Vn-i) = v. (2) 

(In quantum mechanics, h (> 0) is Planck's constant, and the eigenvalue 
parameter —X = Eis the energy.) 

Semiclassical methods are small-^ treatments or expansions for quantum 
problems such as eq.(|l]). Generically, they deliver factorially divergent se- 
ries in h, which at best qualify as asymptotic approximations for the true 
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quantum solutions. Nevertheless, a pioneering work of Balian-Bloch |I| has 
established that fully exact semiclassical formalisms are also achievable in 
principle. This idea was then concretely implemented into an exact WKB 
method for ID problems like eq.(|l|), using resurgence analysis and (general- 
ized) Borel resummation of ^-expansions in this approach, however, 
a gap still remains in the proof of a required resurgence theorem (||^, thm 
1.2.1 and its Comment), and moreover, some of the most explicit calculation 



steps were reached only for very special cases 0-||TT . 

The subject of this review is a newer exact semiclassical (WKB) method 
12[-|T^ more directly based upon functional relations over the original pa- 



rameter space, in the line of Sibuya's framework [0. This treatment is easier 
to formulate at once in broad generality — still within ordinary (linear) dif- 
ferential equations (i.e., ID quantum problems). Moreover, this approach re- 
veals striking unexplained coincidences (currently centered on homogeneous 
differential equations) with Bethe-Ansatz methods for some classes of exactly 
solvable 2D problems: certain statistical-mechanical lattice models and con- 



formal quantum field theories ||16||-p^ 



As far as we know, exact semiclassical formalisms can be modeled on 
old-fashioned (asymptotic) ones. Our own "cooking recipe", in fact, is to 
"exactify" (rewrite in exact form) a conventional semiclassical framework: 
here, a WKB formalism. We will further specialize to a polarized Sturm- 
Liouville problem, given by eq.(|l]) on a half-line only: [0, +oo) without loss 
of generality, and with a choice of Neumann or Dirichlet boundary conditions 
at g = 0. Accordingly, we will recall some basic WKB results (Sec.l), then 
the exact ingredients that we incorporate (Sec. 2), and thereafter how they 
together combine into an exact semiclassical formalism, first for eigenvalue 
problems (Sec. 3), and finally for eq.(|l|) in full generality (Sec. 4). We intend to 
stress the governing ideas throughout, at the expense of technicalities which 
received emphasis in our previous works. 



1 Standard asymptotic inputs 

Our essential starting semiclassical ingredient is asymptotic WKB theory in 
one dimension: specifically, solutions of eq.(|I]) ("wave functions", in quantum 
language) admit the WKB approximations 

V'±(g) ~nA(g)-'/'e±H'nA(</')d5'^ ^here (3) 
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nA(?) = \ly{.<l) + A (= classical forbidden-region momentum). (4) 

Moreover, as is well known, eq.(|^) holds in a compound asymptotics sense: 
not only for — 0, but also, when the potential is polynomial, at fixed h for 
|nA(Q')| — * oo, which includes |A| oo or |g| oo, at will. Henceforth, we 
choose to scale out h to unity. Then, two distinct asymptotic consequences 
can be drawn from eq.(|^). 

1.1 Wave- function asymptotics {q +oo) 

We first introduce a "classical" q +oo expansion, with the notation 

[r(g) + A]-+V2 .;,A)g'^-^^ (a = |, ^ - 1, ^ - 2, . . .); (5) 

we will mostly need the leading (3a- with a > — 1 (which are independent of 
A, with one singular exception: when = 2 !) and at s = 0, which we then 
denote Pa-{v). 

We now return to the quantum problem: following Sibuya []r3|, a reces- 
sive or subdominant solution ip\{q) of eq.([l|) on the half-line can be fully 
specified by a decaying WKB behavior like (|^) in the q +oo direction. 
The reexpansion of that WKB expression in q using eq.(|^) yields 




only a constant normalization factor C remains to be assigned, (e.g., C = 1 
in ll^). Then, the solution ipx gets uniquely fixed through the condition (H). 

1.2 Eigenvalue asymptotics (A —oo) 

Here we consider the potential V^d^'l), thus extended to the whole real line 
by mirror symmetry. It is classically confining, and gives rise to a discrete 
quantum energy spectrum S (for the variable —A). Let the eigenvalues be 
{Ek}k=Q,i,2,... in increasing order: then, Ek t +oo, and the parity of the index 
k moreover matches the parity of the eigenfunctions; equivalently, even (resp. 
odd) k label the Neumann spectrum (resp. Dirichlet spectrum S~) on 
the half-line [0, +0x3). 
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Now, the WKB approximation (^) for the eigenfunctions entails a semi- 
classical formula for E = E^, called Bohr-Sommerfeld quantization condi- 
tion: ^ ^ 

^ ^ k + - for integer k — > +00, (7) 



l{p2+V(\q\)=E} 27r 2 

valid asymptotically for large k, i.e., E — > +00. If that condition is further 
expanded to all orders in k~^, then reorganized in descending powers of E, 
it takes the form [Q: E = E^ obeys 

1 12 

K{v)E'^ ^^^2 integer k +00 (a = /x, l^- j^i ^~ n'"'^' 

(8) 

where the leading exponent is /i = |+;^ (with = (27r)^^ /|p2_,_|g|jv=£;} pdg, 
the classical action for the homogeneous V case). Notes: the ba{v) are poly- 
nomial in the {vj}j<{^-a)N, and actually parity-dependent as well (unless V 
is an even polynomial) but only from a = —3/2 downwards; proper quantum 
corrections (beyond eq.(|^) only contribute to the orders a < — /i, which are 
not critically needed here. 



2 Exact inputs 

We now explain how the preceding WKB approach can be made exact by 
incorporating a single exact ingredient, essentially: that the Wronskian of 
two solutions of eq.([|) is a constant! (precisely, it will be the computable, 
nonzero Wronskian of two particular solutions) . 



2.1 Spectral determinants 



The main exact spectral function we will seek is the canonical or zeta- 
regularized determinant of the eigenvalues -D(A), formally Ylf,{X + Ef^). It can 
be rigorously defined in terms of a spectral Hurwitz-zeta function, J2k{^ + 
Ek)~^, as an entire function of order /i in the variable A (and implicitly entire 
in v), through [|| 



Yi>^ + Eky 



(9) 



We will actually need the fixed-parity determinants -D^(A) (for the sepa- 
rate spectra S^), plus exact limit formulae (of Euler-Maclaurin type) which 
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effectively rebuild each determinant from its spectrum, as 



logD±(A) = lim I J2 log(^fc + A) + ^ \og{EK + A) (10) 



{a>0} 



log ^if - - 



Here the second line expresses counterterms, which have a semiclassical char- 
acter: each of them corresponds to a term of the diverging (as E — > oo) initial 
part of the expansion (^, so as to achieve a finite K —>■ +00 limit. Eq.(p!OD 
amounts to limit-product representations of D^, which are rigidly deter- 
mined by tfie zeros alone (unlike Hadamard products, which retain extra 
free factors). 

2.2 The basic exact identities 

By suitably integrating the Schrodinger equation along the half-line, we can 
obtain a pair of identities relating the spectral determinants to a subdominant 
solution i'xi^q): 

D-{X)^MO), Z^+(A) = -V^UO), (11) 
subject to this explicit specification for the normalizing factor C in eq.(^): 

\ogC^l\-i2log2)p.,{s;v) + dJ^-^^)] (12) 



1 - 2s 



s=0 



Eqs. (pTIJT^) answer a "central connection problem", as they provide some 
kind of data at g = for a solution that was fixed by its q —>■ +00 behavior. 

Remarks: in [jl2|, |13|, we erroneously presumed that C always equals unity, 
whereas this is only guaranteed when /3_i(s, tT) = (which happens not so 
rarely, however); the general case was fixed in (Corrigendum) and |14 



In the latter reference, the result (|T^ was tied to the fact that "classical 
determinants" -D^(A) can also be defined, again by eq. ([TT|) but replacing the 
exact ipx by its WKB approximation (from eq.(^). 

2.3 The conjugate problems 

As with roots of an algebraic equation (which are better studied all together), 
we will not treat the Schrodinger equation (|l]) in isolation, but together with 
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a properly defined set of "conjugate" equations, involving the continuation 
of eq.([l|) to the whole complex domain, and analytic dilations controlled by 
an angle ^ = An/^N + 2). 

The Grst conjugate Schrodinger equation is eq.(0) with potential V^^^ and 
spectral variable A^^l, defined as |]15| 

e-'^ y(e-i^/2 g), aW = e"'^ A. (13) 

The iterated rotations (by £ip) likewise generate a sequence of conjugate 
equations, which repeats itself with a Gnite period L given by 

+ 2 generically 

N (14) 
— + 1 for even polynomials V{q). 

The distinct conjugate equations are thus labeled by [i] = i modL, and 
likewise for all related quantities: e.g., the (A^ — l)-uplet of coefficients v 
from eq.(0) induces a sequence of conjugates -ut^l 



2.4 The Wronskian identity 

Now comes the key exact input, derivable from all preceding considerations. 

Let ipx\(l) be the recessive (for q +oo) solution identically defined for 
the first conjugate equation; then, by counter-rotation (analytic dilation), 

^a(5') == "^aw (s"^''^ I?) is a solution to the original equation (|I]) decaying in 
an adjacent Stokes sector relative to the first solution ip\{q). Here, no Stokes 
phenomenon has to be resolved: both solutions can be asymptotically tracked 
in a common direction, e.g., q +oo, and their (constant) Wronskian ex- 
tracted in this limit |]T3[, giving 

^A(g)^l(g) - mq)Mq) = 2C^'^Cie'^/^e'^^~^^'^y' . (15) 

If the basic identities ([TT|) are substituted into (|15|), with C^^l = 1/C following 
from eq.(0), this Wronskian formula also reads as 

e+W4 /}+(e-i^ A; v^^^)D-{X; v)-e-''^/^ D+(A; v)D-{e-''^ A; t;™) = 2ie'^^-^("^)/' 

(16) 

(reinstating the parametric dependence of the determinants upon v). 

In standard approaches, this explicit Wronskian is viewed as an elemen- 
tary but auxiliary input, and the main dynamical information is sought in 
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Stokes multipliers, namely Wronskians of nonadjacent solution pairs. For in- 
stance, a typical "basic" Stokes multiplier (one involving a next-to- adjacent 
solution pair) has a formula (unnormalized) analogous to eq.(|l^), 



(17) 

However, such multipliers are accessible only through connection formulae 
involving Stokes phenomena, hence in general they are not explicitly known. 

By contrast, our claim is that complete exact information can be squeezed 
out of the elementary Wronskian (|15]) or (|T6D alone. 



3 Exact analysis of eigenvalue problem 

We now argue the sufficiency of eq.(|T6|) as exact input to recover the spec- 
tral determinants (or the eigenvalues of V^(|g|)). Remark: such an exact 
WKB analysis for the Airy functions, using the potential |g| was a 

basic milestone to understand the general case: its validity testified that a 
kink singularity is harmless in the potential (although it looks pathological 
from a naive complex- WKB viewpoint, i.e., worse than any complex-analytic 
singularity) . 

At first sight, the single functional equation (0) seems definitely incom- 
plete to determine the two unknown functions -D^(A; v). Formally, its general 
solution should allow one function to remain arbitrary. Concretely, we may 
even write the full system of all L distinct conjugates of eq. (p!6D , then solve 
it either for in terms of D~ , or vice- versa, but in any case one function 
out of the pair will stay wholly unconstrained. 

Remark: invoking Stokes multipliers is of no help in this perspective; for 
instance, eq.(p!7D adds one equation indeed, but also one unknown function, 
C(A; v) itself. 



3.1 An exact quantization condition 

First of all, we shift to the spectra as basic unknowns, given that the de- 
terminants are fully recoverable from them through eq.(^). For instance, 
let us seek the odd spectrum £_ (the even spectrum can be separately pro- 
cessed likewise). Setting A = —Ek with k odd already reduces eq.(|l6D to a 
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(2/71) ArgD 3 _ 




2 3 4 
(2/7t) log ID'I 



ArgD 



Figure 1: The exact quantization formula (20) for the odd eigenvalues of the 
homogeneous quartic potential, V{q) = q^: the curve {D~ {q™/^ E\^)} !^e>q} inter- 
sects the straight line {ArgI?~(A;0) = 27r/3 modvr} precisely at = E}., with 
k odd (complex-logarithmic coordinates are used; the -E-scale along the curve is 
labeled in italics). 



multiplicative form: 

- e-'^/^ D+{-Ek]v)D- {- e-'^ Ek]^^^) = 2ie'^'^-^^'^^/'^ {k = 1,3,5,...)- 

(18) 

Then, dividing this by its first conjugate formula, we manage to eliminate 
the even spectrum altogether, obtaining 

Z)-(_e+i^E,;^hi]) - ^ (fc- 1,3,5,...). (19) 

Finally, by taking logarithms (fixed by analytical continuation from the ori- 
gin in {E , v)-space) , we end up with a formula fulfilled by the eigenvalues, 
which has an outer Bohr-Sommerfeld form (namely, F{E) = k + 1/2) but is 
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completely exact: for k odd, E = obeys 

^ArgD-(- e-^ E- ^+'^) - ^(5^,{v) = k + \- (fc = 1, 3, 5, . . .)• 

(20) 

This is an exact spectral formula for indeed... but in terms of (the 
determinant corresponding to) another similar spectrum, the first conjugate 
£^3. So, this step still ends with two unknowns (spectra) for a single con- 
straint (the quantization condition (pO|)). 

Special example (homogeneous case): if V{q) = then /5_i(0;0, A) = 
Sn,2^/'2, all conjugate spectra coincide, and the left-hand side of eq.(pOD 
reduces to (2/7r)ArgL)-(- e"'"^ 0) + 6^,2 E/2; fig.l illustrates this in the 
case N = 4. 



3.2 Complete set of exact quantization conditions 

Next, we observe that each of the conjugate potentials V^^^{q) for given i 
can be handled wholly in parallel to the original potential. Consequently, its 
spectra S± also obey exact quantization conditions extending eq.(pOD and 



its even-fc analog. They are now complex in general, and involve the two 
adjacent spectra (of the same parity) S±^^^ and S± ^\ as 



logZ^±(-e-'^E;#+i]) -logD±(-e+'^E;t;[^-i])l - {-lY^(3_,{v) 



to which must be appended all the conjugates of eq. ([TO|) , 



logD±(A; t;M) = hm ( log(4'' + A) + J \og{E^^ + A) (22) 



{a>0} 



log^'-- 



a 



The latter serve to eliminate the determinants from eqs.(0); moreover, 
eqs.(^) imply extra conditions through the coefficients ba{v^^^), namely the 
asymptotic (Bohr-Sommerfeld) formulae extending eq.(|]) to all conjugate 
spectra, 

ba{v^^^)[El!Y -k + - for integer k ^ +oo : (23) 



these have to be obeyed, simply for the definiteness of eqs.(^). But the 
present exact quantization conditions (|21|) are far too imphcit to include the 
semiclassical ones, eqs.(^), as a limit: like the Wronskian identity (^6|) (of 
which they are just an offspring), eqs. (|2T|) explicitly incorporate only the 
degree N plus one invariant, P-i{v), out of the whole dynamics (they are 
"quasi-universal"); whereas the totality of the coefficients Vj of the potential 
are reflected in eqs.(p3D (we expand these at least over {a > — /i}), and 
nowhere else. So, the asymptotic Bohr-Sommerfeld formulae ([23| ) need to be 
asserted independently here; they now come as crucial boundary conditions 
(for k = +00) to eqs.(pH). 

Remark: for the difference term in brackets (left-hand side of eg . (|2l|) ) , 
eqs.(^) reduce to a convergent series (except for = 1 or 2), namely 
J2k[^og{E]^^ — e~"^ E) — \og{E]^^ — e"'"'''' E)]. The counterterms are then optional 
(but still very helpful to accelerate the convergence of this series). 

Now, because the conjugates form a finite cyclic set, each full system 
([21|^) (for a fixed parity) precisely gives L exact equations for L unknown 
spectra. Equivalently, we have reached a potentially complete and closed 
system for each global unknown or S'_, which we define as: £^ =^ U[^] ({^} x 
giV£;M^^ the labeled union of all conjugate spectra (rotated by e'^*^ for sheer 
convenience) . 

Just as the homogeneous-potential cases correspond to 2D integrable 
models and their Bethe-Ansatz exact "solutions" [0-[|TB], the exact quan- 
tization conditions (21-p^) can also be viewed as Bethe-Ansatz equations 
for the general polynomial- potential case. In the same line, we may then 
say that the Schrodinger eigenvalue problem (|l]) gets "integrated" by that 
Bethe-Ansatz system, presuming that the latter can be solved uniquely and 
effectively. 

The challenge posed by the eigenvalue problem in eq.(|I|) hence gets dis- 
placed onto the new system (pi|-pB|), which is of a very different nature: it is 
nonlinear, selfconsistent (or bootstrapping), and above all, purely discrete, no 
longer differential at all — albeit still (countably) infinite-dimensional. Such 
elaborate novel features make us unable to further discuss the solvability of 
that system on rigorous terms (even for the simplest cases of homogeneous 
potentials, e.g., V{q) = q^); e.g., we cannot answer basic questions such as: 

- can one compute backwards: do the exact quantization conditions (plD 
imply the Wronskian identities ((p!6|) and its conjugates) for the determinants 
reconstructed by eq.(p2D, and perhaps even the eigenvalue property ([|) ? 
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- can the system (PT|-P^ have a (reasonably?) unique solution; optimally, 
can it be turned into a fixed-point problem for some contractive mapping(s), 
in the form M^{S^} = 

3.3 Numerical tests 



We are nevertheless able to tackle the resolution of eqs.(|2l|-|23D numerically. 
This system can be approximated on computer by i^- dimensional mapping 
schemes (with K large). Our implementation exploits a "division of labor" 
which appears to exist within that system (though not as sharply as we now 
state it): for k — > +00 {k > K numerically), where asymptotic semiclas- 

sical theory is accurate enough to be used alone, the eigenvalues i?|f' get 
essentially determined by the explicit asymptotic Bohr-Sommerfeld formu- 
lae (|23|), and this enacts the boundary conditions; the finite-fc (or k < K) 
eigenvalues then adjust themselves selfconsistently in response to the implicit 
exact quantization conditions (|2l|) , to which they are more sensitive. 

We have thus tested homogeneous potentials (up to q^^'^) on the one hand 



11], and moderately inhomogeneous quartic and sextic potentials on the 



other hand |]12|. The results strongly support the most favorable picture 
proposed above, namely that nonlinear dynamics of eqs.(^T]-23) can be con- 
tractive. Under straightforward iterations, our numerical mappings indeed 
displayed geometric convergence, and towards the correct quantum spectra 
(up to finite- ii' effects). Within each fixed degree N, the homogeneous case 
looked the most contractive. For quartic potentials V{q) = q'^ + vq"^, we 
have validated such convergence in the range — 10 < f < +5: as we de- 
creased V away from (double-well regime), we quit at t> = —10 after seeing 
nothing worse than a graceful numerical degradation (a gradual decrease of 
final accuracy); by contrast, in the opposite direction (single- well regime), 
iteration schemes could easily be made stable for moderate v, but each of 
them (so far) rather abruptly "derailed" into numerical instability at some 
point, V ^ +5 being our highest result. So, we cannot yet see the expected 
smooth transition towards the harmonic regime as f — > +00 within those 
calculations. 

3.4 Analytical generalizations 

The numerical difficulty just mentioned led us to further study the case 
= + vq^ purely analytically. The v 00 regime is singular indeed 
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as it induces a discontinuity (a jump from 4 to 2) in the degree N, the 
controlhng factor of the exact analysis. But far from reveahng any intrinsic 
singularities, the exact WKB formalism proved regular also in the large-u 
regime (IQ, Sec. 3). For instance, it supphed the large-f behavior for the 
(zero-energy) spectral determinants Qi^(v) *= det^(— d^/dg^ + + vq'^): 



Qi^iv) ~ e-^'''/3 det±(-dVdg2 + vq') [= e-^^'^r^ ^;±V«^^^__:j ^ (34) 



plus many exact results for these functions Qi''^(f) (mostly around f = 0, 
whereas standard perturbation theory operates around v = +00). 

Extension is straightforward to general binomial potentials, i.e., to the 
functions det=^(-dVdg^ + q^ + vq^) (M < A^). The t;-roots of those de- 
terminants solve the generalized eigenvalue problems 

-riq) + [q'' + vq^'Mq) = 0, (25) 

in which the coupling v is now taken as spectral parameter whilst the original 
one is frozen (A = 0), and eq.(^) is posed over {q > 0} just as before, 
with Neumann or Dirichlet conditions at g = 0. The cases where N = 
2M + 2 are special; for even M and selected v, these give supersymmetric 
potentials. Exact WKB analysis still works here (fl^) Sec. 4; |18[): it yields 



exact quantization conditions for v, which degenerate to a totally explicit 
form when N = 2M + 2 (compare with eq.(plD for ^ = 0): 



--p-i{v) = - -^T^v = k + -± ;' for k = "'^f -. (26) 



TT 



1 , A^-2 



0,2,4,. 



N + 2 2 2(A^ + 2) ^'3,5,. 



The usual harmonic oscillator spectral problem, formerly an isolated ex- 
ception, is now seen as the first in an infinite sequence of solvable cases 
(iV = 2,6,10,...). 

A wholly novel exactly solvable family is also uncovered by that exact 
WKB analysis, namely the generalized spectrum of eq. (|25|) over the whole 
real line when N = 2M + 2 with M odd. (For even M, this spectrum just 
consists of both parity components of eq.(p6D.) For odd M, the full spectral 



determinant over the whole real line comes out as |jT4| 
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implying the exact quantization condition 



(2M + 4)(n + 1/2) (nGZ) 



if M is odd. 



(28) 



Here, even the lowest case M = 1 does not reduce to a classic solvable case. 



The corresponding eigenfunctions have also been recently investigated [|T9|. 
Furthermore, both degenerate quantization conditions i 



2^) get identi- 



cally reproduced by the leading semiclassical Bohr-Sommerfeld formula (|^: 
i.e., as for the harmonic oscillator, "semiclassical quantization is exact" (but 
for those generalized eigenvalues only). 



4 Exact wave- function analysis 

4.1 A solution algorithm for the Schrodinger equation 

We finally show how the exact treatment can solve the differential equation 
(d) itself, i.e., compute the recessive solution ^ipx^q) for instance, in just a 
slight extension of the previous exact eigenvalue calculations: it suffices to 
use the basic identities (pAj) in reverse. The choice of origin being immaterial, 
those can be written over the half-line [q, +oo) as 

Mq) = D;iX + Viq)), i^'M = -^,+ (A + V{q)) (29) 

where -D^(A) are the spectral determinants of the potential V^(-) defined by 
\V{-) — V{q)\ over the half-line [g, +cxd), and by mirror symmetry on the 
complementary half-line. 

Thus, for every fixed g, ip\{q) can be computed from the Dirichlet spec- 
trum £q^^ of the potential V^(-) (and likewise, ipx from the Neumann spec- 
trum). That spectrum in turn obeys an exact quantization formula, now 
g-dependent, hence can derive from an explicit fixed-point equation of the 
form A4~{S' _} = Once this spectrum has been attained, possibly by 
iteration, one extra application of the reconstruction formula like eq. ( |T0| ) for 
the determinant D~ can supply any value of this function, in particular the 
one required by eq. (p^) . 

We emphasize that all distinct points q get handled completely rigidly 
and autonomously here, invoking no propagation from one point to the next. 
In this sense, the whole procedure that we have described "integrates" (one 
solution for) the Schrodinger equation (|l|). 
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Figure 2: Exact- WKB calculations for a ground-state eigenfunction tp{q), in the 
potential V{q) = — 5g^. Dashed curves (using the right-hand vertical scale): 
the potential V{q) and its ground-state energy level Eq ~ —3.41014, to show the 
classically allowed {{V{q) < Eq}) and forbidden {{V{q) > Eq}) portions of the q- 
axis; vertical lines mark the two inner turning points. Continuous curve (using 
the left-hand vertical scale): eigenfunction ^l^{q) resulting from a computer inte- 
gration of eq.(|l|) (using the NAG routine D02KEF), and normalized according to 
eqs.(I^Jl^) with C = 1 here (which is not normalization!). + : exact- WKB values 
ip{q), reached by iterative calculations of the determinants D~ of eq.(^), inde- 
pendently for each sampled g- value (numerical instability disrupts this calculation 
scheme beyond \q\ 1.7). 

4.2 Numerical tests 

We have validated this exact WKB algorithm numerically upon several exam- 
ples with quartic potentials, and for solutions ipil) in fully quantum regimes; 
we provide two illustrations. 

Fig. 2 displays a ground-state wave function (for a symmetric double- 
well potential, featuring = 0). Now, inasmuch as the present approach 
works for all recessive {q — * -|-oo) solutions alike, it also cannot directly 
detect which ones will end up being square-integrable for g — — oo as well; 
the easiest way to pin down an eigenfunction is then to preassign —X = E 
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Figure 3: Exact- WKB calculations for a non-L^ recessive solution ipo{q), in the 
potential V{q) = — 2q and at energy E = 0. Continuous curve (now using 
a logarithmic vertical scale): analytical solution ipoiq), normalized according to 
eqs.(I^JT^) with C = \/2 here. + : exact- WKB computed values ipo{q)^ as in fig.2. 



itself as the eigenvalue. So, the ground-state eigenvalue Eq —3.41014) is 
determined separately first, e.g., by means of the basic exact quantization 
conditions (|2T|) ; the resulting number A = —Eq is then fed into the calculation 
of every value ipx^q) according to eqs. (pUj) and the associated g-dependent 
quantization conditions. As in the eigenvalue problem, a numerical instability 
abruptly sets in outside a parameter range which is {|g| < 1.7} here (but 
the contracting factors of our iterations are numerically ~ 0.67 for most of 
the computed g-points, and do not exceed 0.86 at the farthest computable 
ones). Furthermore, the stability range encompasses two turning points {q ~ 
±0.90267), which clearly do not affect the exact WKB calculations. 

Fig. 3 shows a contrasting example: a non- square- integr able solution, now 
for an asymmetrical quartic potential featuring a nonzero /9-i(s; v) (= 2s — 1). 
The chosen problem is eq.(P5|) with = 4, M = 1, and with v = —2 not 
belonging to the generalized spectrum (given by eq.(^)): in this problem, 
a supersymmetry (on {q > 0} with Neumann boundary conditions) makes 
the fully normalized solution available in closed form. The iterations again 
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display contraction factors ~ 0.67 (at all plotted points). 



4.3 Concluding remarks 

We first restate our two major open issues (which might be fruitfully merged). 
On the one hand, the nonlinear infinite-dimensional structure of eqs. (|2T|^3D 
is entirely uncontrolled at a rigorous level. On the other hand, the analytical 
coincidences observed with exactly solvable 2D theories |jl6|]-|l^ are in utter 



want of fundamental explanations (and, hopefully, of generalizations too). 

We are also confident that such approaches will extend beyond polynomial 
potentials. On the basis of existing results {^^, resp. pO|)' believe that 
exact-WKB resolution methods should next be sought for rational, resp. 
trigonometric-polynomial, potentials; a later realistic and desirable target 
could be all Heun-class equations. 
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